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ABSTRACT 

During rapid solidification processing (RSP) the amount of liquid undercooling Is an important 
factor In determining microstructural development by controlling phase selection during nucleatlon and 
morphological evolution during crystal growth. While undercooling Is an inherent feature of many 
techniques of RSP, the deepest undercoolings and most controlled studies have been possible in care- 
fully prepared fine droplet samples. From past work and recent advances In studies of nucleatlon 
kinetics It has become clear that the Initiation of crystallization during RSP is governed usually 
by heterogeneous sites located at surfaces. With known nucleant sites, it has been possible to 
Identify specific pathways of metastable phase formation and microstructural development In alloys. 
These advances have allowed for a clearer assessment of the interplay between undercooling, cooling 
rate and particle size statistics In structure formation. New approaches to the examination of growth 
processes have been developed to follow the thermal behavior and morphology in small samples In the 
period of rapid crystallization and recalescence. Based upon the new experimental information from 
these studies, useful models can be developed for the overall solidification process to Include nu- 
cleation behavior, thermodynamic constraints, thermal history, growth kinetics, solute redistribution 
and resulting structures. From the refinement of knowledge concerning the underlying factors that 
govern RSP a basis Is emerging for an effective alloy design and processing strategy. 


Introduction 


The main attention in rapid solidification (RSP) methods has been focused upon the attainment of 
high cooling rates in the range of 10^ to 10° °C/sec. With these methods small liquid volumes are 
quenched either in contact with a substrate as in melt spinning and directed energy beam processing or 
in an inert gas as in atomization. However, the most important impact of RSP can be related directly 
to the solidification structures and the new oppor-^unities for structure control by innovative alloy 
design and processing strategies. The wide range o, reported microstructural variations encompass not 
only equilibrium phase mixtures with refined microstructural scale, but also novel microstructures 
such as supersaturated solid solutions, metastable intermediate phases and amorphous solids. 

The existence and understanding of novel microstructures representing non-equilibrium phases 
requires a consideration of the level of undercooling at the onset of nucleatlon in addition to rapid 
growth kinetics. In fact, the available free energy for non-equilibrium phase formation is directly 
related to the amount of undercooling. With increasing levels of undercooling the product phase 
selection involved in nucleatlon can be expanded to involve an increased variety and to allow competi- 
tive growth kinetics to play an important role in the evolution of different microstructural morpholo- 
gies. Accordingly, the examination of solidification in undercooled liquids can contribute to the 
basic understanding necessary for an effective microstructure control and alloy design in the ootimiza- 
tion of RSP. 

Experimental Techniques for High Undercooling 

Even with liquid metals of high purity the attainment of a high level of undercooling is precluded 
usually by the cat£' "tic action of residual impurities or container walls. In order to observe exten- 
sive undercooling i s necessary to circumvent the catalytic effects of the active nucleants. 

An effective experime. al approach that may be applied to yield large undercooling prior to solidifica- 
tion involves the slow cooling (10-30°C/min) from the melt of a dispersion of ctabilized fine metal 
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droplets [1], By starting with high purity metal and dispersing the liquid Into a large number of 
small droplets (<20y) within a suitable medium, the catalytic effects of active nucleants may be res- 
tricted to a small fraction of the total droplet population. A suitable means of carrying out the 
nucleant isolation operation is illustrated schematically in fig. 1. 

In considering the droplet approach to undercooling, a distinction should be made bet'-een the 
several modifications of this technique which have been employed by a number of investigators in the 
past. For exoi.ile in the droplet substrate method [2], a small sample of liquid metal is placed on an 
inert glass substrate in a chamber with ah inert or reducing atmosphere. With this approach the onset 
of solidification is determined by visual observation of the change in surface reflectivity of a metal 
droplet which is about 50-100p in diameter, In the entrained droplet technique [3] a bulk alloy 
sample is equilibrated in a liquid-solid two phase field with the temperature and alloy composition 
selected so that only a small fraction of the volume will be liquid in contact with a primary solid 
solution phase, as a potential heterogeneous nucleant. In addition, several other procedures have 
been developed for generating liquid metal droplet dispersions which exhibit a propensity for substan- 
tial undercooling. These droplet forming techniques involve exploding wires [4J, shotting [5] and the 
containerless solidification in drop tubes [6]. While the majority of droplet undercooling schemes 
have a provision for some method of temperature measurement, the undercooling tendency of liquid drop- 
lets in several cases including all atomization methods is judged usually from microstructural observa- 
tions. 

With bulk samples, the removal of active nucleation sites usually involves a physical separation 
of the sites from the melt or a chemical treatment either to incorporate or to deactivate possible 
catalytic sites in the presence of an inorganic containment layer. These approaches which include a 
bulk sample encased in organic glasses [7] as well as a bulk melt levitated by induction [8] indicate 
that some remarkably large levels of undercooling (i.e. 0.2 Tm) may be achieved, The inorganic glass 
encasing is believed to promote undercooling not only by isolating the liquid from possible catalytic 
sites on container walls, but also by a possible scavenging action on nucleating centers that may be 
distributed within the volume of the melt. It i- worthwhile to note that in principle the undercooling 
attained in properly conditioned bulk samples should approach the levels in droplet samples. Besides 
these specific undercooling techniques, it is clear that atomization of liquids is a commercial RSP 
method which incorporates an effective nucleant isolation operation as well as rapid quenching. 




Nuclaation Kinetics In Undercooled Liquids 


In the present discussion the main Interest 1n nucleatlon kinetics will be focused on the possible 
relationships between undercooling and cooling rate during the solidification of a distribution of 
liquid droplets. These relationships may be examined by nucleatlon rate calculations, but the computed 
results are highly dependent on the specific values assigned to the kinetic rate parameters. Reliable 
Information Is available for some parameters, however, 1n most cases the kinetic parameters are assign- 
ed values based upon theoretical models which require experimental verification. To minimize the 
uncertainty associated with calculations based model dependent parameter values. It is of interest to 
examine the range of droplet nucleatlon kinetics In a general manner. With this approach It is possible 
to deduce several useful guiding relationships for rapid solidification behavior. 

Nucleatlon of a crystalline phase In an undercooled liquid pha.‘-e may be classified Into two mech- 
anistic types. When nucleatlon occurs without an association with any catalytic nucleatlon sites. It 
car. be classified as homogeneous nucleatlon. In the vast majority of experimental conditions nuclea- 
tlon Is Initiated at some catalytic nucleatlon site as heterogeneous nucleatlon. For a heterogeneous 
nucleatlon, an expression for the nucleatlon rate, J, has been developed on the basis of homogeneous 
nucleatlon theory [9] as: 


J = Hg exp C-AG(n*)/kT] (1) 


where S2a proportional to the number of liquid atoms in contact with unit area of the catalytic 
surface as well as the jumping frequency of atoms, AG(n*) is the activation barrier for nucleatlon and 
kT has the usual meaning. For nuclei having the form of a spherical sector, 

AG(n') = (16/3) f(e)/(AG^)^ (2) 

where a is the Interfacial energy between the nucleus and the liquid, AGy is the driving free energy 
for nucleatlon of unit volume of the product phase and f(0) = (2-3 cos 6 + cos^ 0)/4 Is a contact 
angle factor. In the descrip<"'on of nucleatlon kinetics during rapid liquid quenching, steady state 
nucleatlon conditions may be considered in general, except for glass forming alloys [9]. 

Nucleatlon Rate Parameters 


The nucleatlon rate Is a relati'»ely steep function of temperature with a steady state magnitude 
determined principally by the value of the exponential terms in AG(n*) at the nucleatlon temperature 
Tn and to a lesser amount by the prefactor. In evaluating the temperature dependence of J usually a 
constant value based on theory is taken for Qa sven though the catalytic site may vary. Similarly, 
little information is available to judge the catalytic potency, f(0) so that a range of values is 
used in calculations. However, the most Important parameters in determining J and hence the under- 
cooling level are AGy and o which have received continued experimental and theoretical study. 

In pure metals, the recent measurements of the heat capacity for undercooled liquids [10] have 
allowed for a direct evaluation of AGy. In general, the experimental values are represented to within 
a few percent by the relation proposed by Turnbull [11,12], 


AG^ = -AK^ (AT/TJ = - AH^ (1-T^) (3) 


where AHy is the heat of fusion per unit volume of product phase, AT is the undercooling and T^ is the 
melting point of the product phase. For alloys the value of AG^ is determined from the composition of 
the nucleating phase at T^. The composition dependence of AGy can be expressed as 

AG^ = [X^ Au^ + (1-X^) Au^]/(? (4) 

A B 

where Ap and Ap are the chemical potential differences between the undercooled liquid and crystal and 
V is the molar volume for the crystalline phase at composition Xft. It is apparent from eq. 4 that AGy 
is maximized when Ap^ = Ap^ so that AGy = Ap^/v. This condition has been proposed by Hillert [13] to 
assess the nucleus composition. However, the favored nucleus composition should yield a minimum value 
of AG(n*). Only if v and a are not composition dependent can the maximum AGy value occur at the same 



composition as that for a minimum In AG(n*). Appa'-.^ntly this behavior may be approached In some alloys. 
For example, recent calculations [14] of the composition dependence of the nucleation temperature, 

T„(x) have reproduced the observed experimental trends [1] which reflect the composition dependence of 
tne llquidus, but the generality of the calculated behavior requires further examination. 

In the study of nucleation one of the most Important parameters is the solid-liquid Interfacial 
energy. Since the nucleation rate is proportional to exp -(o)3 as shown in eqs. 1 and 2, it is very 
sensitive to even slight errors in the o value [15]. A direct measurement of a is very difficult 
even with a large uncertainties because of sensitivity to impurity effects. Usually values of o are 
calculated from classical rucleatlon theory by assuming that experimentally obtained maximum undercool- 
ing values are associated with homogeneous nucleation. Unless the maximum undercooling value is 
proven to refer to homogeneous nucleation, the estimated a value is expected to underestimate the 
true value. There have been a number of theoretical approaches to study the structure and properties 
of the solid-liquid interface on the basis of various model constructions. The models range from the 
limiting cases in which o values are calculated from enthalpic bond breaking arguments [lo] or derived 
solely from a basis of entropic considerations [17] to intermediate approaches where a weighted con- 
tribution of enthalpic and entropic terms [18-20] is Incorporated into the evaluation. At the present 
time, the theoretical values appear to underestimate cr compared to values calculated from nucleation 
undercooling results [15]. Certainly, an improvement in the Independent experimental determination 
of o would be invaluable in allowing for a clear evaluation of theoretical models and more reliable 
nucleation calculations. 

Effect of Cooling Rate on Nucleation Kinetics 

Since the onset of nucleation in a highly undercooled liquid is expected to be a sharp function of 
temperature it is useful to consider a limiting case for the effect of cooling rate on nucleation 
temperature. In general, for a given nucleation kinetics the critical condition to observe nucleation 
experimentally in a time, t, in a droplet of volume, v, and catalytic surface area, a, is given by 
Jvt 1 for homogeneous nucleation and by Jat ~ 1 for heterogeneous nucleation. From eq. 1, the 
heterogeneous nucleation time can be represented as a function of temperature as; 

In t = -ln[a (n^)]+ 16 TO3f(e)/[3 k (AH^)^ T^ 0>T^)^ T^] (4) 


where a r,g represents the heterogeneous nucleation site population associated with a droplet of volume 

V. 


Furthermore, the rate dependence of the nucleation temperature which represents the shape of the 
T-T-T diagram associatid with a specific nucleation kinetics is represented in general as 

•jTTiffr ‘ „3 LfJTTn J 


where the possible timperature dependence of ilg, AHy, a and f(6) has been neglected in order to ex- 
amine a limiting value and minimize the uncertainty of model dependent values. For homogeneous 
nucleation f(S) = 1, but for a heterogeneous nucleation Q<f(e)<1. Therefore, T-T-T diagrams for 
heterogeneous nucleation can have different breadths representing a different temperature dependence 
of the nucleation rate for different heterogeneous sites. In eq. 5, dTr/d(1nt) becomes infinite at 
Tr = 1/3 as the nucleation time is minimized. Therefore, T^ = 1/3 is believed to estimate a lower 
bound to the onset of sensible nucleation of a crystalline phase during slow cooling [21]. If the 
undercooling limit of the droplet sample represents the solidification of the nucleant-free droplet 
population, the crystallization temperature is determined either by the catalytic nature of the drop- 
let surface coating or by the intervention of homogeneous nucleation kinetics. Also, during cooling 
the operation of a given catalyst kinetics and the kinetic competition between different catalysts and 
different product phase structures will be influenced strongly by the active catalyst distribution at 
the undercooling limit. In order to reach this limiting undercooling the size refinement of droplets 
becomes one of the most important factors. Further analysis on the distribution of active nucleants 
among droplets then becomes necessary to understand the undercooling behavior of droplet samples. 

Nucleation Kinetics and Size Distribution Statistics 


As long as a single nucleation kinetics with a steep temperature dependence is operating the 
effect of droplet size on the nucleation temperature is relatively small. In fact, if the size 
distribution is the only factor associated with the nucleation temperature range during cooling, a 



sharp nucleat^on peak over a narrow temperature range (<10°C) 1s expected even with a sample of a 
rather broad droplet .:l 2 e distribution. However, droplet emulsion samples of pure metals and alloys 
often show a broad crystallization peak over a wide temperature range (>20°C). In this cane the 
operation of a single nucleatlon kinetics Is unlikely but the broad crystallization peak can be 
associated wit!, multiple nucleatlon events continuously activated during the cooling process [21]. 
Morever, at a gl/en cooling rate, nucleatlon of the crystalline phasse In each droplet Is associated 
with the catalytic potency of the most active nucleant. As a result, nucleatlon at different levels 
of undercooling due to crystallization over a wide temperature range can result In the formation of 
different structures In different droplets. Even If the same structure Is nucleated at different 
levels of undercooling, this structure will not necessarily be retained completely In the final solid- 
ification pattern. For example, the final product structure can develop variations due to a different 
thermal history associated with the recalescence period following nucleatlon at different undercooling 
levels In separate droplets. 

An example of the nucleatlon of different phases due to crystallization behavior over a broad 
undercooling range may be Illustrated by considering tie careful studies on Fe-Ni alloys reported by 
Cech [22]. In fine droplet samples, (<40p), Cech demonstrated the presence of droplets of a metastable 
bcc structure in addition to droplets of the equilibrium fee structure. Size distribution analysis 
results showed that .small droplets less than 13u are more likely nucleated at high undercooling since 
the major fraction had a metastable bcc structure. S1r>ce there were also droplets of the metastable 
structure within the large size grouping (d > 20y ), It Is necessary to consider the mode of nucleant 
distribution among droplets. When all possible nucleatlon sites which can be activated a lower levels 
of undercooling than that necessary to retain the metastable phase are considered to be nucleants, a 
Poisson distribution function may be introduced to describe the random distribution of these nucleants 
among droplets [21,23]. In this case the fraction of droplets free from active nucleants Is represented 
by 


X = exp [- m (v)] (6) 

where m(v) is the average number of active nucleants in each droplet of volume V. Depending on the 
distribution mechanism of nucleants during the atomization process m(v) Is proportional either to the 
volume or to the corresponding surface area of the droplet. Since a nucleant distribution proportional 
to the droplet volume would result In a rapid decrease in the yield of metastable product with Increas- 
ing droplet size. It appears that for the Fe-NI case, a surface area dependent nucleant distribution is 
most probable. A further application of the Poisson distribution has been reported for the solidifica- 
tion of Pd-Si droplet samples [24]. In this case the droplet fraction which solidified continuously 
as glass represents X In eq. 6 which can also be written in terms of droplet diameter, d, as X = exp 
[-(d/do)’’]. When m(v) is proportional to the droplet surface area, n=2 for a heterogeneous surface 
nucleatlon. For heterogeneous volume nucleatlon n=3 and for homogeneous nucleatlon n=4.6. As in other 
droplet studies the sol irificatlon behavior of Pd-Si droplets was described most closely by a hetero- 
geneous surface nucleatlon kinetics. 

Metastable Phase Formation 


Thermodynamic Features 

A large number of examples of metastable phase formation have been demonstrated during rapid 
quenching experiments [25]. While the possible reaction paths leading to metastable phase formation 
in alloy systems are numerous a consideration of thermodynamic constraints can be useful in providing 
guidelines and limiting bounds on the reaction sequences in undercooled liquids. For example the free 
energy relationships that are relevant to the generation of solid solutions are illustrated in fig. 

2. Two features are of interest in this case. For the formation of an equilibrium solid solution of 
composition Cq at temperature Tfj an overall free energy change of aGi is Involved In reaction. However, 
in the early nucleatlon stage of reaction when a small amount of a forms from the liquid the driving 
free energy is maximized for a with a composition C] rather Cg. When the molar volume of solution does 
not change appreciably with composition, the C] composition may be determined by constructing a tangent 
to the a free energy curve which is parallel to a tangent to the liquid free energy. As noted previou- 
sly, this consideration indicates that a with solute content different from Cq can develop initially 
during nucleatlon [13]. Another feature of interest in fig. 2 concerns the bound set by the Tg dashed 
curve between the liquidus and solidus [26], which represents the locus of equal free energy between 
li,uid and solid. Formation of a product with the same composition as the liquid is excluded at 
compositions beyond the Tg limit. 

A thermodynamic consideration is of most value when the analysis also can provide an estimate of 
potential metastable intermediate phases that can enter in the kinetic competition. For example, in 
alloys the requirement that the free energy, A3, exhibit an overall decrease with reaction can lead 
to several interesting possibilities for phase reaction sequences [26] as illustrated by one example in 
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F1g. 2 Free energy relationship between under 
cooled liquid and a phase. 



Fig, 3 Free energy relationship between under- 
cooled liquid and several crystalline 
phases [26]. 
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Fig. 4 Schematic illustration of a possible reaction path for primary phase catalysis 
of a metastable solid solution (a) and a metastable y intermediate phase (b). 
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fig. 3. In this case the equnibrlum solidification of liquid of composition Cq would yield a two- 
phase mixture of a with composition C^, and 6 with composition Cg. However, if the nucleation kinetics 
are favorable, an intermediate reaction step can occur involving the formation of the metastable y 
phase with composition Cq anc free energy reduction AG(L-»y. )• The presence of y Phase affects tiie 
reaction sequence in that a phase cannot form since the free energy change ) is positive. 

Only after the y phase transforms into S phase can the a phase form with a free energy reduction. 

When there are several metastable phases possible, the number of intermediate reaction steps can In- 
crease and present interesting opportunities for structure control during both solidification and sub- 
sequent solid state processing. One way of summarizing the potential for structure modification 1n 
alloys is to construct a metastable phase diagram [25,27]. 

Competitive Kinetics Analaysis 

Since the driving force for the formation of a metastable ptiase is related directly to the amount 
of undercooling below the melting temperature of the phase, a high level of undercooling can expand the 
product selection involved In nucleation to an Increased variety of types of structures. This is 
illustrated by the nucleation transition from fee to bcc In Fe-Ni alloys [22], stainless steels [28] 
and tool steels [29]. Other examples have been reported in Pb-Bi, Bi-Cd and Pb-Sn alloys as well as 
Ga and B1 [1,10,21]. For Pb-Bi alloys a multiple phase selection transition has been observed In 
droplets [1]. As the level of undercooling increases an equilibrium eutectic mixture Is replaced by 
a supersaturated (hep) phase which in turn is superceded by a metastable Intermediate phase at the 
highest undercooling. In addition for Pb-Sn, Sn-Ge and Cd-Sb alloys the operation of the nucleation 
product transition has been related to the catalytic effect of specific heterogeneous sites [21,30]. 

One type of effective nucleation site Is a primary solid solution which can participate in catalysis 
reactions such as those Illustrated in fig. 4 [21]. When the a phase is a potent catalyst, a continued 

growth of a 1s possible during cooling to form a supersaturated a product as In fig. 4a. When a is 

strongly catalytic in promoting the nucleation of a metastable intermediate y phase as Indicated in 
fig. 4b, the formation of an equilibrium S phase may be bypassed as the liquid Is undercooled below 
the metastable y phase melting point. Thus, for an alloy that may exhibit metastable products, it 
Is possible to select a particular phase by controlling the degree of undercooling prior to solidifica- 
tion and by providing a high density of favorable catalytic sites. 

In order to understand the kinetic competition between different phases duri»"g nucleation In an 
unuercooled liquid, the comparison of thermodynamic stability may not be sufficient. In terms of 
overall nucleation kinetics, the driving free energy is not the only parameter of importance, but also 
nucleation site density and the magnitude of AG(n*) associated with a given phase are critical. As an 
illustration of the application of nucleation kinetics to reveal the relationships between cooling rate, 
undercooling and phase selection, it is useful to consider the competition between two nucleation 
kinetics during continuous cooling of a droplet in several different situations that are ^llustrated 
in fig. 5. As an initial example, when two catalysts of different nucleation site density and potency 
are in the same droplet (8 and C), the undercooling limit of the droplet at Trl is deteimined by 
catalyst C at cooling rate t] . At cooling rate t2, however, the undercooling of the droplet increases 
to Tf2 by circumventing the catalytic effect of C. Thus at a continually increasing level of cooling 
rate, there can be a significant Improvement of undercooling as the catalytic effects of highly potent 
nucleants are circumvented. During this process the controlling nucleation kinetics involves catalysts 
of decreasing potency. When two different phases are competing with each other on the same catalytic 
surface present in a highly undercooled liquid, according to nucleation theory, the magnitude of the 
activation free energy barrier becomes the most dominant factor. Therefore, in order to favor nuclea- 
tion of a metastable phase at a given cooling rate the activation energy barrier for the metastable 
phase should be lower than that for the equilibrium phase. In terms of eq. 1 this condition may be 
expressed as 


0, M f(e') .,3^/0, 

LM ^ < 0 f(6) 

(1-T;)^ (AH^)2 (1-T^)2 

where aim and 9' are the interfacial energy and the contact angle between the metastable phaseand the 
catalytic surface, AHy'" is the heat of fusion per unit volume of the metastable phase and T^' is the 
relative nucleation temperature of the metastable phase. It has been shown in eq. 5 that the breadth 
of the T-T-T diagram associated with different nucleation kinetics is proportional to (AHv)‘^Tm/o‘^T(0) 
When the melting temperatures of two different phases nucleated by the same catalysts are very close 
and the metastable phase is dominant, the condition given by eq. 7 may be approximated by 






Fig, 5 T-T-T diagram representing different nucleation kinetics that may 
occur during a continuous cooling in undercooled droplets. 
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In this case, the breadth of the T-T-T diagram for the nucleation of the metastabla phase (A) is larger 
than that for the equilibrium phase (B), but the noses of these two diagrams are placed at the same 
position as illustrated in fig. 5 schematically (A and B), Unless there are nucleants of higher 
catalytic potency for nucleation of the equilibrium phase (for example C in fig. 5), the equilibrium 
phase is always masked by nucleation of the metastable phase regardless of the cooling rate. In fact, 
the droplet emulsion technique which disperses potent nucleants of small population into a minor frac- 
tion of the droplet sample can be viewed as eliminating the possible activation of nucleation kinetics 
C in the major fraction of the sample. Alternatively, if the activation barrier for A is higher than 
that for B, the equilibrium phase will nucleate at all cooling rates. In addition, it is also possible 
for the nucleation of the equilibrium phase to be favored by an increase of the cooling rate if the 
nucleation site density for the equilibrium phase (curvw D) is greater than that for the metastable 
phase. For these conditions an increase in cooling rate can result in a decreasing yield of a metast- 
able product. In fact, a transition in nucleation from surface dependent to volume dependent kinetics 
with increasing cooling rate such as C to D in fig, 5 can yield a large increase in rate and may con- 
tribute to a propensity for multiple nucleation as reported in quenched Al-Si alloy droplets [31]. 
However, other factors involving the specific thermal history of droplets are also likely to contribute 
to the development of multiple nucleation. 

In summary, the examination of the simple competitive kinetics examples reveals several types of 
general behavior. First, for a single operating kinetics an increase in cooling rate translates 
usually into some increase in undercooling. With several competing kinetics the increase in under- 
cooling with cooling rate may not be uniform. For several competing kinetics and possible product 
structures, an increase in cooling rate can allow for nucleation of a metastable phase to be favored 
or for an increase in the yield of a metastable phase product, but this is not a general rule. Depend- 
ing on the density of catalytic sites, it is possible for an increase in cooling rate to result in no 
change in phase selection or to result in a decreasing yield of a given metastable phase in favor of 
another metastable product or even an equilibrium phase product. From these examples illustrating the 
possible synergistic effects between high undercooling and rapid quenching to yield unique product 
structures, it is clear that attention to catalyst potency and site distribution as affected by surface 
conditions during atomization can be as important as cooling rate in optimizing RSP treatments. 




Heat Flow and Pure Metal Solidification 


In solidification of small droplets, two extreme heat flow situations may arise. In the first 
one, the solidification front moves Into a liquid at Its melting point or above, and the latent 
heat of fusion must be withdrawn through the solid. The second situation Involves solidification 
Into a supercooled melt. Nucieatlon occurs only after substantial supercooling, and Initial 
solidification then occurs primarily by flow of heat Into the liquid Itself, therby resulting In 
droplet recalescence. 

For solidification without undercooling, the Newtonian approximation Is applicable when the 
dimensionless Blot number, B1 , Is lower than 0.1, where B1 = h ro/k^. The Blot number Is the 
ratio of the surface conductance or heat transfer coefficient, h to the metal conductance{k^/rQ) : 

Is the thermal conductivity of the metal and rg Is the droplet radius. The heat transfer 
coefficient Is determined from the Nusselt number [32]: 


Nu ■= 2 h r^/k^ -2 + 0.6 Re'^^Fr’^^ (9) 

where Re and Pc are Reynolds and Prandtl numbers respectively and k^ Is the thennal conductivity 
of the fluid. In the case of solidification of metal droplets In a quiescent liquid or gas, 
the Nusselt number reduces to 2, and the Blot number 1s simplified to k^/k . In the case of pure 
tin droplets emulsified In oil (polyphenyl ether), the Bloc number Is about 4.4 x 10"\ since 
k/! Is 0.112 kcal/m hr ®C and k^ Is 25.3 kcal/m hr ®C. Thus, the dominant resistance to heat flow 
In this example occurs at the surface and no significant temperature gradients can develop inside 
the droplet In absence of undercooling. Similarly In most If not all gas atomization processes, 
heat transfer Is limited by the solid-gas heat transfer resistance. 

Thus, the simple Newtonian model Is widely and correctly used for single phase cooling and 
solidification without undercooling [33]. Table 1 shows cooling rates for aluminum gas atomized 
droplets calculated In this way. Table 2 shows comparable cooling rates for atomized or 
emulsified metals In several liquid media. 

When supercooling occurs In atomized or emulsified droplets, it Is not obvious that the 
droplet will have a uniform temperature distribution during the recalescing period of solidification; 

Table 1. Calculated Cooling Rates of Pure Aluminum at 660'’C During Gas 
Atomization (In Helium at 300®K); (®K/sec) 


Relative 

Gas Velocity 

Particle Diameter (vm) 

mach 

V„(ni s"M 

10 

100 

0 

0 

7.21 X 10- 

7.21 X 10^ 

0.1 

97 

1.00 ;< 10^ 

1 .56 X 10® 

0.3 

290 

1.22 X 10' 

2.26 X 10® 

0.5 

■ 483 

1.36 X 10^ 

2.69 X 10® 

1.0 

965 

1.56 X 10^ 

3.53 X 10® 


Table 2. Calculated Cooling Rates of Pure Metal Droplets (lOpm dia.) 
Cooled In Different Media at Low Reynolds Numbers 


Metal and Quench Medium 

Cooling Rate at Equilibrium 
Meltinq Point (°K/sec) 

Aluminum in Helium gas at 27°C 

7.21 X 10® 

Aluminum in Argon gas at 27°C 

9.68 X 10® 

Aluminum in Liquid Gallium at 30°C 

2.12 X 10® 

Tin in Argon at 27°C 

4.57 X 10® 

Tin in Oil at 100°C 

2.28 X 10® 

Tin in 011/CCl4(l/2) at -25’’C 

3.77 X 10® 

Iron in Argon at 27°C 

1 .15 X lO® 

Iron In Glass at 1000°C 

3.99 X lU^ 

Iron in Molten Salt at 482°C 

3.11 X 10® 



I.e., the temperature of the growing Interface may rise faster than the average temperature of 
the droplet. In this case, the validity of the Newtonian approximation during recalescence will 
depend not only on the Blot number, but also on Interface velocity and Interface morphology. 


Levi and Mehrablan [34] have considered In some detail the solidification of undercooled 
droplets of pure metals, assuming a planar (ron»dendr1t1c) Interface and the Interface kinetic 
relationship suggested by Turnbull [35] and modified by Cahn [36] for the steady state continuous 
growth of a planar front. Figure 6 Is an example of reff ''ts of their computer calculations for an 
aluminum droplet with a single nucleatlon event. For conditions chosen, significant temperature 
differences develop within the liquid during recalescence. Of course, as pointed out by Levi and 
Mehrablan, the magnitude and longevity of Internal temperature differences depend sensitively on 
the number of nuclei growing, and on whether or not the Interface Is dendritic. In the case of 
multiple nucleatlon, the scale of the heat transfer problem Is reduced from that of the droplet 
radius, rg, to that of half of the Interparticle spacing dn)x/2. In the limit of lacy dendrites 
(or, dendrites that become unstable and break Into Individual particles), the heat transfer scale 
Is similarly reduced. Figure 7, also from Levi and Mehrablan, shows schematically how decreasing 
d^j^ althers the "enthalpy path" during solidification. Decreasing dm^ moves paths more toward the 
top and right of the diagram and here recalescence becomes more rapid, more adiabatic, and more 
Newtonian. 



Figure 6. Thermal profiles for a recalescing 
pure aluminum droplet, lOpm d1a. 
Advancing liquid solid front shown in 
insert which denotes fraction solid 
and interface veloc1ty(cm s‘M> f'"cm' 
Levi and Mehrablan [34]. 
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Figure 7. Effect of the ratio K/Bi on the 
thermal history of undercooled 
droplets for a fixed nucleatlon 
temperature. Note that for typical 
atomization conditions (k/B1>100) the 
external cooling is not able to slow 
down recalescence to a significant 
extent, from Levi and Mehrablan [34]. 


We do not yet have experimental measurements of growth velocities or Interface morphology In 
undercooled pure metal droplets. We can, however, turn to other studies where growth of crystals 
Into pure undercooled liquids has been measured, and where this growth is apparently dendritic. 
Perhaps the highest measured crystal growth velocities are these of nickel and cobalt Into their 
pure undercooled melts, reported by Walker [37], Colligan and Bayles [38]. In their experiments, 
columns of the bulk melt were undercooled by varying amounts up to 250-300°C and then seeded at one 
end. The velocity was determined from the times at which the freely moving dendrite tip passed 
thermal sensors sited along the column. Similar experiments have been performed by Glicksman et al 
[39] on sued nonitrile, and successfully analyzed by combining heat flow and interface stability 
theory. By addition of interface kinetics to this analysis, Coriell and Turnbull [40] showed 
excellent agreement of theory with experiment, for the undercooled nickel experiments described 
above. These experimental studies, as well as interface stability theory suggest solidification of 
undercooled pure metals (and alloys) should be dendritic. Nonetheless, there are Indications that 
undercooled alloy droplets at high undercooling solidify in a non-dendritic fashion. Perhaps this 
is because of transient effects, or perhaps because the wavwlength of the stability is long compared 
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with the droplet size [41]. 

Alloy Solidification 

In an earlier section, the possibility has been described of achieving non-equilibrium phases 
through undercooling (e.g. obtaining a BCC phase In a normally austenitic steel). And of course, 
glassy structures can also be obtained. In this section, we focus on structure and segregation 
In alloys v^here the expected phases are formed (I.e. those obtained In usual solidification process). 

A number of Investigators have shown that solidification in bulk Iron, nickel, and copper 
alloys Is dendritic at moderate undercoolings, with the dendrite structure becoming gradually a 
simple rod-like structure with Increasing undercooling. Then, at a critical u.i jprcool 1 ng which 
often seems to be about IOO-200'’C grain size drops abruptly; an example Is shown In Figure 8 [42]. 

For undercoolings greater than the critical, the dendrite structure Is completely absent, replaced 
by more-or-less spherical ''elements'' which are the separate grains [7,42]. Both the dendrite arm 
spacing In moderately undercooled specimens, and the size of the spherical elements are strongly 
dependent on cooling rate after recalescence. Hence, it Is clear that the structures we see are 
strongly affected by ripening and are not therefore the actual structures that grew during recales- 
cence. The abrupt grain size reduction has been althernately explained by an acoustic nucleatlon 
mechanism, and by a ripening or dendrite remalting mechanism. 

Especially at Intermediate to rapid rates of heat extraction during and/or Immediately after 
recalescence, a number of Investigators have observed solute rich cores within the dendrites or 
"elements" In undercooled alloys. Figure 9 Is an example from work of Wu [43] showing high solute 
cores In 5mm dia. Fe-262N1 alloy that was undercooled 90°C. The solute rich cores certainly 
formed during recalescence. Perhaps they grew In diffusionless fashion below the Tg, but the 
probable remelting, coarsening, and solid state diffusion occuHng during and after recalescence 
did not completely eliminate these traces of th'i mglnal structure. 

Related structures have been obtained by Levi and Mehrabian In aluminum rapidly solidified 
droplets. Homogeneous solute rich material Is found In a portion (or all) of the droplets when 
a given droplet is not completely homogeneous jc contains also a solute depleted zone and dendritic 
or cellular region as well. The homogeneous solute rich region must have formed at some temperature 
below the thermodynamic Tg temperature, and grew In an apparently non-dendritic manner, at some 
critical temperature during recalescence, solute began to be rejected, eventually rendering It 
unstable and resulting In the cellular structure at the upper right. 
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Figure 8. Grain size versus undercooling In 
a bulk undercooled melt, from 
Kattamis and Flemings [42]. 



Figure 9. Solute rich cores In undercooled 
Fe-25%N1 alloy, AT=90°C, 700X, 
from Wu [43]. 


With sufficient undercooling and sufficiently rapid heat extraction, supersaturated structures 
and non-equilibrium phases can be retained to room temperature. Kelly, Vander Sande and Cohen [28] 
have obtained a homogeneous (non-cel lular), non equilibrium BCC structure 1n atomized 303 stainless 
steel powder in diameter up to about 70 microns. Figure 10 shows similar homogeneous structures In 
316 stainless steel obtained by Macisaac [44] 1n thesis work at MIT. These latter structures are 
of both the equilibrium FCC phase and the non-equilibrium BCC phase; they were obtained in an 
emulsification of steel in glass. These structures also appear to have grown In non-dendritic manner. 
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Figure 10. (•) Optical Micrograph of crost*sact1ontd and atchad 316 BCC and FCC particles 

after undercooling. Fee particles are light In color, (b) Backscattered 
electron Micrograph of the same particles, from Nacisaac [44]. 


The attributes of high undercooling and rapid heat extraction are optimized In emulsions 
when the emulsion Itself Is rapidly cooled by external means. One way to do this Is shown In 
Figure 11 [45], where an emulsion of a metal In a carrier fluid Is Injected Into a cooler miscible 
fluid that Is being vloorously agitated. This Is a technique developed by M.G. Chu In his doctoral 
thesis work at MIT. Figure 12 shows two structures of Sn>5%Pb alloy Illustrating the Importance of 
this external cooling. When the alloy Is emulsified In such a way that It can achieve an under* 
cooling of 100*C or more at slow cooling rates, and then Is rapidly cooled using the apparatus of 
Figure 11, a fully homogeneous, non-dendr1t1c structure Is obtained, at least at 1000X, Figure 12a. 
This Is In spite of the fact that the maximum equilibrium solid solubility of lead in kin Is 2.5t. 
When, however, these samples are more slowly cooled, or reheated, or held for an extended period of 
time at room temperature, substantial precipitation of lead rich phase occurs. Figure 12b. In the 
case of Chu's work, the effect of the quench was probably primarily In preventing solid state 
decomposition. In other cases (e.g.) quenching liquid Iron In molten salt. Table 2, heat extraction 
rates might conceivably be Increased enough by the cold salt to Influence solidification kinetics. 
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Figure 11. Diagram of apparatus for rapid 
quenching of molten emulsified 
droplets, from Chu et a1 [45]. 
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Figure 12. Microstructures of emulsified 

droplets of Sn-5XPb: (a) Rapidly 
cooled to -25*C, (b) Rapidly cooled 

and then held one month at room 
temperature lOOOX. 
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Quito complex solidification structures can be obtained in undercooled alloys when more than 
one phase forms during recalescence. Chu has studied a wide variety of such structures in tin-lead 
alloys; the structures obtained depend on alloy, undercooling and rate of heat extraction during 
and after recalescence. with the particular emulsification technique he used, and a moderate 
cooling rate {10°C/min), he found the nucleation temperatures plotted in the phase diagram in 
Figure 13a. A tin rich phase forms at Tn and a lead rich phase at T,!,. Porepezko [21] previously 
had found similar nucleation temperatures. Figure 13b, and by using a different emulsification 
medium, he was able to get the lower nucleation temperatures also shown in Figure 13b. 

Solidification of hyper-eutectic tin-lead alloys was found to occur in the following sequence: 
(a) nucleation at T/, of a lead rich phase, (b) dendritic growth of this phase, first during recales- 
cence and then during subsequent cooling Tn, where a tin rich phase nucleates, (c) growth of the tin 
rich phase during rP''»Mscence and during subsequent cooling. Note that the recalescence and solute 
distribution during g.owth of the second phase can be expected to influence the first phase to form. 
In addition the intervention of the formation of a metastable phase prior to the tin rich phase due 
to catalysis by the lead rich phase [16] can also influence the overall solidification sequence and 
microstructure. Figure 14 shows examples of solidification structures of such an alloy. Figure 14a 
is a strucutre that was undercooled and rapidly quenched. Lead rich dendrites (white) first formed 
and tin rich phase (dark) later filled in vdthin the interstices. Both the tin rich and lead rich 
phases shown substantial second phase with themselves, p*'esumab1y from solid state precipitation. 
Figure 14b is a sample that was allowed to nucleate on slow cooling and then rapidly solidified. 

The dendritic structure formed as before, and the final structure is very little different. The 
sample allowed to cool slowly (lO®C/m1n) however, shows massive ripening of the original lead rich 
phase. Figure 14c. 

Three different tin-lead structures of eutectic composition are shown in Figure 15. The first 
is comparable to the slowly cooled Sn-45%Pb alloy discussed above, in that it shows primary, ripened 
lead rich phase (with entrapped and precipitated tin rich phase). Figure 15a. The second shows the 
same alloy rapidly cooled. Figure 15b. Apparently the rapid cooling suppressed the nucleation 
temperature of the lead rich phase, T„ below that of the tin rich phase, Tn. The third structure. 
Figure 15c, from work of Cooper, Anderson and Perepezko [46] is similar to 15b, but somewhat coarser. 
In this case the nucleation temperature T^ was depressed belo'^ Tp, by a modified emulsion; this 
depression is shown in figure 13. The second two eutectic st uctures are similar to those obtained 
in conventional rapid solidification processes, e.g. wheel casting [47]. 


Recalescpo '.t «' - surements in Droplet Solidification 


£w_ . [4bJ and Wu [43] 1" 
seperatf! theses at MIT are attempting 
to gain some experimental under- 
standing of the "solidification 
paths" of undercooled alloys by 
looking at thermal behavior and 
structure of 5mm dia. metal droplets 
undercooled in molten glass. Their 
results to date illustrate well 
how much we do not understand yet 
about this process. 

As one example, the rate of 
recalescence of the droplet surface 
increases dramatically with in- 
creasing undercooling, Figiure 16. 

We believe that this cannot be a 
heat flow caused phenomenon since 
experiments on pure nickel show 
orders of magnitude faster recales- 
cence. Thus, a preliminary con- 
clusion Is that these rates of 
recalescence (at least at tempera- 
tures approaching and above the 
calculated Tq) are limited by 
solute transport. The Tq tempera- 
ture in Figure 16, incidentally, 
was calculated from a regular 
solution model . 


Compostli'jn of trad (4l*7*Pb) Compoaition of laad (al%Fb) 
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Figure 13. Nucleation temperature in supercooled Sn-Pb 
alloys. T_ is nucleation temperature of 
tin rich pnase and T^ of lead rich phase. 
Diagram on left is from Chu et a1 [45] and 
on the right is from Cooper et a1 [46]. 

Exact amounts of undercooling depend on 
emulsification procedure. The two different 
set^ of curves in diagram on the right result 
from the different emulsification procedures. 







Flqure 14. Microstructurts of Sn*451Pb alloy dropletk, 7000X. (a) Rapidly qutnchtd 

(b) Slowly cooltd to Just balow Tp and quanchad (c) Slowly coolad to room 
tamparatura. 
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Figure 15. Microstructures of enulsifled droplets of eutectic Sn-Pb alloy, (a) Slowly 
cooled, moderate undercooling, (b) Undercooled and rapidly solidified, 

(c) Slowly cooled, high undercooling. 


The photomicrographs in Figure 17 
provide some Insight Into factors In- 
fluencing this recalescence, by showing 
the very different structur*-*, st different 
undercoolings. At mode.'ate urdercool Ing 
where temperature Is estimated to be 
above the thermodynamic temperature, 
recalescence Is very slow, and a two 
phv'se structure result: that has a 
coarse dendritic appearance. Figure 17a. 
Initial results show the core of this 
"ciandrlte" Is cf appreximately Cq 
composition, the periphery Is depleted, 
and the Interdendr It^c regions appear 
to be somewhat above C^, possibly It 
grew as a single phase^lnltially 
Into solute redistribution at least 
during later stages of growth. If 
so, we would expect re-melting on 
a fine scale within the dendrites 
as It recalesces above some tempera- 
ture approaching T ; the liquid 
pools remelting could lead to the 
finely distributed second phase 
within the "dendrite". 



Figure 16. Recalescence curves from 5mm 

droplets of Nl-Sn eutectic alloy 
at several undercoolings, from 
Ewasko [48]. 


The sample at the Intennedlate undercooling In Figure 17 show a fine grain structure with a 
multiplicity of growth sites, each of which has a central region (two phase) that may well have 




InItUlly grown In • di ffui1onl#»$ iwnnor t>«1ow T_ *nd 1«t#r partially remclted during rtcalescenca. 
Th# fine scale of this structure due to the multiplicity of growth sites would explain the more 
rapid recalescence, even above T . This qualitative description of solidification behavior for 
such an alloy Is largely conjecture, but the approach suggests an experimental direction toward 
an eventual fuller thermal and solidification description of undercooled alloys. 



Undercooling 70*C Undercooling 120*C Undercooling 380*C 


Figure 17. Structure of 5nwn droplets of Nl-Sn eutectic alloy soHdlfed at 3 different 
undercoolings, from Wu [43]. 

Conclusions 


In examining the nucleatlon and solidification of undercooled liquid metals It Is apparent that 
the experimental and theoretical understanding has made good progress since the Initial pioneering 
work of Turnbull. From the present experience with droplets studies It has been possible to Identify 
some general guidelines for the Interplay between process variables such as cooling rate, under- 
cooling, alloy composition and particle size and the nucleatlon phase selection Involved In RSP. 

With droplet samples and controlled undercooling conditions new understanding has been gained Into 
the development of solidification microstructure and the Influence of thermal history on the final 
microstructural morphology. This area of research In the study of RSP, we believe, has a great 
potential for the further development of nucleatlon and solidification theory. It also has the 
potential for the development of new processes of engineering Importance. In this regard It seems 
possible that we can someday learn to enhance undercooling with controlled surface treatments In 
existing RSP methods, and perhaps, as well, find ways to bring other processes such as emulsification 
to the point of engineering usefulness. 
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